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I.  INTRODUCTION 


In  developing  computer  programs  to  numerically  solve  the  Navier-Stokes 
equations,  the  purpose  of  the  computation  must  be  clearly  kept  in  mind.  In 
the  Air  Force  Flight  Dynamics  Laboratory,  our  purpose  is  to  provide  design 
information  on  "non-linear"  aerodynamic  phenomenon  for  aircraft  that  per¬ 
form  throughout  the  flight  corridor.  This  translates  into  requirements  for 
a  computer  program  which  can  solve  the  time  averaged  compressible  Navier- 
Stokes  equations  (with  a  turbulence  model)  in  three  dimensions  for  general¬ 
ized  geometries.  The  intended  application  of  the  results  then  establishes 
the  priorities  in  addressing  critical  issues. 

In  our  investigations  of  viscous  flows,  several  problem  areas  keep 
recurring.  They  are  as  follows: 

1)  Grid  generation  for  arbitrary  geometry 

2)  Numerical  difficulties 

3)  Turbulence  models 

4)  Accuracy  and  efficiency 

5)  Smearing  of  discontinuities 

A.  GRID  GENERATION  FOR  ARBITRARY  GEOMETRY 

It  is  generally  accepted  that  viscous  flow  problems  require  a  surface- 
oriented  coordinate  system.  Also  for  arbitrary  geometries,  automation  of 
a  numerical  transformation  (as  opposed  to  an  analytic  transformation)  is 
necessary.  In  addition,  some  optimization  of  the  distribution  of  grid 
points  throughout  the  flow  field  is  necessary  to  economically  solve  prac¬ 
tical  problems.  Conceptually,  this  implies  that  higher  order  derivatives 
(in  the  transformed  plane)  of  the  primary  dependent  variable  be  minimized.^ 


The  distribution  of  the  grid  points  greatly  influences  the  requirement  of 
the  number  of  field  points  necessary  to  achieve  a  desired  accuracy.  Con¬ 
siderably  more  attention  is  needed  in  this  area  to  Improve  the  economics 
of  the  viscous  flow  computations. 

B.  NUMERICAL  DIFFICULTIES 

This  is  a  "catch  all"  term  to  cover  the  reasons  a  program  "bombs  out". 
Given  a  proven  algorithm  and  an  experienced  user  with  a  properly  formulated 
problem,  program  failures  are  still  common  during  the  inital  phase  of  the 
investigation.  The  problems  are  most  frequently  due  to  large  truncation 
errors  which  eventually  swamp  the  true  solution.  The  cause  of  the  problem 
is  that  the  grid  cannot  truly  be  established  until  the  flowfield  is  deter¬ 
mined.  A  redistribution  or  increase  in  the  number  of  grid  points  often 
permits  success.  Artfully  changing  the  damping  coefficients  in  the  region 
of  discontinuities  has  also  been  successful.  In  addition,  alternate 
approaches  for  expressing  the  boundary  conditions  can  have  a  dramatic 
effect  on  the  success  or  failure  of  a  problem.  A  requirement  exists  for 
a  method  in  which  the  flowfield  modifies  its  own  numerical  grid  where 
needed  (adaptive  mesh  generation).  Also,  additional  program  guidelines 
are  needed  to  insure  a  more  robust  code. 

C.  TURBULENCE  MODELS 

In  time-averaging  the  Navler-Stokes  equations  Information  is  lost. 
Information  must  be  re-inserted  into  the  governing  equations  by  resorting 
to  experimental  observation.  The  engineer  needs  empirically  determined 
transport  properties  to  proceed  with  the  numerical  computation.  A  large 
body  of  data  exists  for  flat  plate  boundary  layers  and  good  correlations 
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have  evolved  which  generally  permit  calculations  to  be  performed  that  fit 

2 

the  data  to  within  +10%  for  skin  friction  and  boundary  layer  thickness 

(see  Fig.  1).  Unfortunately,  the  agreement  for  the  pressure  gradient 

case  is  not  nearly  as  good.  Higher  order  closure  schemes  have  not  greatly 

improved  the  prediction  capability.  There  is  a  need  for  the  measurement 

of  turbulent  Reynolds  stresses  under  pressure  gradient  for  a  wide  range 

of  flow  conditions  to  permit  correlations  comparable  to  the  flat  plate 

case.  Without  this  data,  progress  in  the  field  will  be  limited. 

Many  skeptics  are  pessimistic  of  our  ability  to  compute  turbulent 

flows  in  the  near  future.  Turbulence  is  felt  to  be  too  complex  and  the 

progress  has  been  slow  in  developing  a  thorough  understanding.  To  counter 

these  skeptics,  an  encouraging  viewpoint  is  offered.  First,  the  good 

design  predictions  of  flat  plate  properties  are  possible  without  fully 

understanding  the  true  mechanism  of  turbulence.  Secondly,  in  some  cases 

it  may  be  possible  to  bracket  the  extremes  of  flows  with  pressure  gradient 

3 

by  computing  the  frozen  and  equilibrium  states  ,  thereby,  providing  useable 

4 

design  information  (Fig.  2).  Thirdly,  remarkable  results  are  possible  in 
the  prediction  of  gross  turbulent  properties  by  simply  treating  the  eddy 
viscosity  Reynolds  number  as  a  constant  (^~^  ■  Ret  *  const).  Turbulence 
is  limited  and  confined,  and  these  approximate  results  are  easy  to  compute; 
the  difficulty  is  in  reducing  the  error  bounds  to  satisfy  the  scientist. 
Fourthly,  in  most  applications,  only  displacement  effects  which  influence 
the  pressure  distribution  (separation  point  location)  are  significant.  Skin 
friction  and  heat  transfer,  which  require  greater  numerical  resolution,  are 
often  of  secondary  importance.  Therefore,  there  should  be  hope  that  it  may 
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be  possible  to  devise  a  turbulence  model  capable  of  producing  useful 
design  information  for  many  engineering  applications. 

One  last  point  concerning  the  future  development  of  turbulence  models. 
The  models  to  date  have  been  analytical  and  thus  ’’analog"  in  nature.  New 
models  have  an  additional  requirement  to  be  compatible  with  numerical  com¬ 
putation.  We  need  something  like  "digital  turbulence." 

D.  ACCURACY  AND  EFFICIENCY 

Accuracy  and  efficiency  should  be  addressed  concurrently  because  of 

their  interrelationship.  Given  a  stable  algorithm,  the  greatest  control 

on  spatial  accuracy  is  the  number  and  distribution  of  grid  points.  Figure 

3  shows  the  error  in  drag  coefficient  vs  number  of  grid  points^  in  one 

coordinate  direction  in  an  airfoil  flowfield.  The  computational  time 
2 

increases  with  N  (for  a  two  dimensional  problem)  and  hence  it  is  very 
expensive  to  obtain  the  last  few  percent  accuracy.  The  accuracy  require¬ 
ments  of  any  design  problem  must  be  very  carefully  defined  in  order  to 
avoid  excessive  computer  cost. 

Once  satisfactory  spatial  accuracy  is  achieved,  a  convergence  criteria 

must  be  selected  which  produces  comparable  accuracy.  A  time  dependent 

approach  is  generally  used  to  solve  the  Navier-Stokes  equations  in  which 

the  computation  proceeds  from  an  arbitrary  initial  condition  until  a  steady 

state  solution  is  achieved.  In  the  past,  several  (maybe  5)  characteristic 

times (||— )  has  been  sufficient  for  the  initial  transient  to  decay.  However, 

based  upon  the  analytical  solution  of  an  impulsively  started  flat  plate, 

—is 

the  error  between  the  transient  value  and  steady  state  decays  as  t  . 

This  slow  convergence  rate  Implies  that  to  cut  the  error  in  half,  the  com¬ 
puter  time  must  be  increased  by  a  factor  of  four^  (for  the  same  At).  (See 


Fig.  4.)  Another  discouraging  aspect  is  that  for  some  flows,  periodic 
values  are  legitimate  steady  state  solutions.  For  example,  subsonic  air¬ 
foils  near  stall  shed  vorticies  in  a  regular  manner^  (Fig.  5).  Computations 
must  be  accomplished  for  many  characteristic  times  to  achieve  mean  and  rms 
values  for  design  application.  Slow  convergence  could  well  be  our  most 
critical  problem  in  our  goal  to  economically  produce  aerodynamic  design 
data. 

Paramount  to  all  of  these  issues  is  the  fact  that  a  good  finite 
difference  algorithm  must  be  used  to  solve  the  governing  equations.  Con¬ 
siderable  success  has  been  achieved  with  MacCormack's  method^  to  solve 
supersonic  viscous  flows.  MacCormack's  explicit  method  possesses  many 
desirable  features  with  the  exception  of  efficiency.  The  CFL  stability 
limit  requires  small  time  steps  where  small  spatial  steps  are  required  to 
resolve  viscous  regions.  To  relieve  this  restriction,  implicit  methods 
have  been  developed  which  are  conceptually  unconditionally  stable.  However, 
our  experience  shows  a  gain  in  efficiency  only  in  the  viscous  region. 
Accuracy  (not  stability)  requirements  in  the  inviscid  region  can  be  achieved 
only  for  the  CFL  time  step.  Hence  the  hybrid  method"*  (explicit  in  the 
inviscid  and  implicit  in  the  viscous  region)  is  at  present  probably  the 
most  efficient  method  available. 

E.  SMEARING  OF  DISCONTINUTIES 

In  examining  viscous  flow  problems,  two  scale  lengths  appear.  One 
is  the  mean  free  path,  X  ~  y-  ;  the  other,  which  is  introduced  through  the 
boundary  conditions,  is  a  characteristic  geometric  length,  L.  One  can 
also  derive  another  scale  length,  5  ~  /LX  ;  which  is  a  combination  of 
the  previously  mentioned  two  lengths.  In  numerically  solving  any  viscous 
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flow  problem,  the  grid  size,  Ay,  should  be  sufficiently  small  to  accurately 

resolve  these  three  scale  lengths  (L,  <5,  A).  This,  of  course,  is  impossible 

to  achieve  in  nearly  any  practical  problem  today.  Slip  lines,  shock  waves 

and  leading  edge  flows  are  examples  where  the  characteristic  lengths  are 

too  small  to  be  honored.  As  a  consequence,  these  discontinuities  are 

incorrectly  computed.  Large  errors  exist  in  the  immediate  vicinity  of 

these  regions  and  numerical  smearing  results.  Based  on  both  wind  tunnel 

and  computational  experience,  it  is  believed  that  these  local  errors  near 

singularities  do  not  totally  invalidate  the  global  results.  Figure  6  shows 

9 

a  Navier-Stokes  computation  of  a  high  speed  inlet  flow  indicating  good 
agreement  with  experiment  with  the  exception  of  the  shock  jump  and  the 
entropy  layer  generated  by  the  cowl  lip  leading  edge.  More  effort  is 
required  to  minimize  the  smearing  of  these  discontinuities. 

F.  APPLICATIONS 

The  principal  application  areas  in  computational  fluid  dynamics  during 
this  reporting  period  can  be  generally  grouped  into  four  categories,  i.e., 
(1)  shock-turbulent  boundary  layer  interactions  (2)  3-D  interactions 
(3)  inlets  and  nozzle  flowflelds  and  (4)  unsteady  viscous  flows.  These 
research  investigations  shall  now  be  discussed. 
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II.  SHOCK  WAVE-TURBULENT  BOUNDARY  LAYER  INTERACTIONS 

The  general  features  of  shock  wave-boundary  layer  Interactions  were 
identified  by  Dr.  Hans  Liepmann  in  1940's.  Analytical  predictions  were 
accomplished  in  the  1960's  by  simultaneously  coupling  the  laminar  boundary 
layer  equations  with  simple  supersonic  wave  theory,  thereby  successfully 
describing  the  interaction  process.  This  elementary  method  encountered 
numerical  difficulties  when  extended  to  turbulent  boundary  layers,  hyper¬ 
sonic  flows  or  complex  geometries.  For  that  reason,  the  time  dependent 

Navier-Stokes  equations  with  an  algebraic  turbulence  model  were  utilized 

XI  12  13  14 

to  solve  the  turbulent  flow  case  '  ’  ’  .  Nearly  instantaneous  success 

resulted  in  that  the  numerical  difficulties  were  overcome  and  the  global 
features  of  the  interaction  described  (Fig  7  and  8).  It  was  discovered, 
however,  that  the  equilibrium  turbulence  model  was  inadequate  to  accurately 
predict  the  size  of  the  separation  bubble  or  the  skin  friction  level  in 
the  reattachment  region.  However,  the  feasibility  to  numerically  compute 
these  flowfields  was  clearly  demonstrated  for  the  first  time.  Current 
research  is  now  being  directed  to  improve  the  computational  efficiency 

and  to  improve  the  accuracy  of  the  turbulence  model  for  adverse  pressure 

..  ►  15,16 

gradients 
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III.  THREE  DIMENSIONAL  INTERACTING  FLOWS 

Early  work  with  the  Navier-Stokes  equations  was  for  simple  2D  geom¬ 
etries.  In  order  to  advance  the  method  for  practical  design  applications 
an  effort  was  undertaken  to  solve  3D  configurations.  Two  approaches  were 
taken.  First,  conical  (quasi-three  dimensional)  flows  were  examined  and 
secondly  the  full  three  dimensional  Navier-Stokes  equations  were  computed. 

Conical  flows  possess  three  velocity  components  as  a  function  of  only 
two  spatial  variables,  and  hence  introduce  some  features  encountered  in 
true  3D  flowfields.  Only  supersonic  inviscid  flows  over  conically  symmetric 
bodies  are  exactly  conical.  However,  in  regions  where  the  radial  derivatives 
of  the  viscous  terms  are  small  compared  with  the  transverse  derivatives  a 

"locally-conical"  approximation  can  be  performed  (Fig.  9).  In  this  manner 

17  18  19  20 

the  flowfields  over  cones  *  and  delta  wings  '  at  angles  of  attack 
were  computed.  Excellent  agreement  was  obtained  for  these  conical  flows 
when  compared  with  experiment  (Fig.  10  and  11). 

Efforts  in  computing  fully  3D  flows  concentrated  on  axial  corner,  shock 
wave-boundary  layer  interactions  primarily  because  of  the  availability  of 
detailed  experimental  flowfield  data.  To  prove  feasibility  comparisons  must 
be  performed  with  experimental  results  which  have  sufficient  detail  to  out¬ 
line  the  essential  flow  features.  The  axial  comer  problem  possessed  such 
data  for  both  laminar  and  turbulent  supersonic  and  hypersonic  conditions. 

Computations  for  these  cases  indeed  proved  the  feasibility  of  the 
21  22  23 

approach  ’  ’  (Fig.  12).  For  the  cases  computed,  the  insensitivity  of 

the  turbulence  model  for  3D  flows  was  observed  (as  opposed  to  the  opposite 
for  2D  flows).  The  cases  investigated  may  not  have  been  severe  enough  to 
display  a  large  discrepancy,  however. 
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The  need  for  improved  computer  efficiency  was  clearly  displayed  in 

computing  these  cases  and  hence  the  research  on  3D  Navier-Stokes  was  di- 

rected  toward  the  use  of  vector  processing  computers.  The  3D  time  dependent 

Navier-Stokes  program  was  converted  to  operate  on  the  CRAY-1  computer  with 

the  assistance  of  the  University  of  Michigan  staff  under  the  direction  of 

Professor  D.  Calahan.  On  September  1978  a  vectorized  version  ran  at  a 
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speed  of  128  times  the  CDC  6600  computer  (Figure  13).  This  accomplishment 
was  a  major  step  towards  the  achievement  of  practical  3D  flowfield  computa¬ 
tions  for  design  application  (Fig  14). 


IV.  INLETS  AND  NOZZLE  FLOWFIELDS 


With  the  foundation  laid  to  compute  viscous  flows  with  shock  waves 
for  simple  shapes,  efforts  focused  on  refining  the  code  to  solve  practical 
configurations  with  complex  geometries.  Computers  possess  insufficient 
storage  to  solve  the  flowfields  for  complete  configurations,  however,  air¬ 
craft  component  flows  can  be  resolved.  Components  requiring  considerable 
design  attention  are  inlets  and  nozzles.  Attention  was,  therefore,  focused 
in  this  area. 

Following  the  guidelines  developed  by  Dr.  Joe  Thompson,  a  numerical 

transformation  from  a  body-oriented  coordinate  system  in  the  physical  plane 

to  a  unit  square  (with  constant  step  size)  in  the  transformed  plane  permits 

25  26  27  28 

a  completely  general  approach  for  resolving  any  configuration  ’  *  ’ 

The  overall  objective  of  this  approach  is  to  establish  a  versatile 

generalized  method,  simplify  the  prescription  of  the  boundary  conditions, 

and  minimize  truncation  error.  The  first  coupling  of  the  compressible 

Navier-Stokes  equations  with  Thompson's  transformation  was  accomplished 

29 

for  a  chemical  laser  diffuser  by  Knight  and  Hankey  .  Others  quickly 

followed  so  that  this  approach  is  presently  the  most  accepted  proce- 
30  31  32  33 

dure  '  ’  ’  (Fig  6) .  In  the  future,  the  method  must  be  refined  to 

include  adaptive  mesh  generation,  local  orthogonality  and  automated  smooth¬ 
ing  of  the  transformation  derivatives. 

General  grid  generation,  stable  shock  capturing  methods,  effective 
turbulence  models  and  stable  boundary  conditions  are  necessary  features 
requiring  attention  in  computing  inlets  and  nozzle  flowfields. 


V.  UNSTEADY  VISCOUS  FLOWS 


Solutions  to  the  Navler-Stokes  equations  are  obtained  by  utilizing 
the  time-dependent  terms  and  computing  from  free  stream  initial  conditions 
until  steady  state  results  are  achieved.  Therefore,  this  approach  has  the 
capability  of  analyzing  some  unsteady  flows  (Fig  15).  A  class  presently 
under  investigation  is  the  self-excited  flow  (i.e.,  buzz,  buffet,  reson¬ 
ance).  These  flows  are  generally  the  result  of  an  unstable  separated  shear 
layer  (Fig  16)  which  grows  until  a  "limit  cycle"  is  approached-1  ’  .  If 

the  reduced  frequency  parameter  is  less  than  unity  it  is  believed  that 
accurate  solutions  are  possible  utilizing  the  Reynolds-averaged  Navier- 

Stokes  equations.  Good  agreement  with  experimental  results  has  been  obtained 

36 

for  a  pressure  oscillation  in  an  open  cavity  (Fig  17).  This  area  shows 
considerable  promise  for  future  payoff  in  analyzing  unsteady  flow  problems 
of  interest  to  the  Air  Force. 


11 


VI.  CONCLUSION 


Although  additional  research  is  required,  we  believe  outstanding 
advances  have  been  accomplished  in  Computational  Aerodynamics.  The  ground 
work  has  been  laid  in  that  stable  algorithms  exist,  mesh  generators  for 
fairly  arbitrary  shapes  have  been  used,  turbulence  models  are  available 
that  produce  global  results,  and  the  use  of  vector  processors  has  demon¬ 
strated  a  rapid  speed  up  of  the  codes.  Although  the  feasibility  has  been 
displayed,  much  improvement  is  needed  in  all  areas  before  the  practically 
of  computational  fluid  dynamics  is  achieved. 
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Figure  1.  Comparison  of  Skin  Friction 
Values  for  a  Flat  Plate 
Boundary  Layer 
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Figure  2.  Comparison  of  Frozen  and 

Equilibrium  Turbulence  Models 
for  Compression  Ramp 


14 


S/L 


GRID  PTS 


o 

oc 

CL 

LU 

Figure  3.  Error  in  Force  Coefficients 

vs.  Number  of  Grid  Points  for 
Inviscid  Flow  Over  an  Airfoil 
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Figure  4.  Convergence  Characteristics  for 
Shock  Wave  Impingement 
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Figure  10.  Velocity  Vector  Flow  Field 
£or  Cone  at  Angle  of  Attack 
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Figure  11.  Cross  Flow  Streamlines  for 

Hypersonic  Flow  Over  Delta 
Wing  a  15° 


Figure  13.  Comparative  Timing  Results 

for  3D  Navier-Stokes  Program 


25 


Figure  14.  Perspective  View  of  Density 
Contours  for  a  Normal  Shock 
in  a  Wind  Tunnel  Diffuser 


Figure  16.  Unstable  Wave  Number  Region 
for  Similar  Separated  Flows 
(Lower  Branch  of  Falkner- 
Skan  Boundary  Layers)  28 
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